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ABSTRACT 

We study the spectral and energetics properties of 47 long-duration gamma-ray bursts 
(GRBs) with known redshift, all of them detected by the Swift satellite. Due to the 
narrow energy range (15-150 keV) of the Swift-B AT detector, the spectral fitting is 
reliable only for fitting models with 2 or 3 parameters. As high uncertainty and cor- 
relation among the errors is expected, a careful analysis of the errors is necessary. 
We fit both the power law (PL, 2 parameters) and cut-off power law (CPL, 3 pa- 
rameters) models to the time-integrated spectra of the 47 bursts, and we present the 
corresponding parameters, their uncertainties, and the correlations among the uncer- 
tainties. The CPL model is reliable only for 29 bursts for which we estimate the vfi, 
peak energy Ep]^. For these GRBs, we calculate the energy fluence and the rest-frame 
isotropic-equivalent radiated energy, iiJ^jso, as well as the propagated uncertainties 
and correlations among them. We explore the distribution of our homogeneous sample 
of GRBs on the rest-frame diagram E'^^^ vs -E-y,iso- We confirm a significant correlation 
between these two quantities (the "Amati" relation) and we verify that, within the un- 
certainty limits, no outliers are present. We also fit the spectra to a Band model with 
the high energy power law index frozen to —2.3, obtaining a rather good agreement 
with the "Amati" relation of non-Swift GRBs. 
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1 INTRODUCTION 

The Swift Mission (|Gehrels et al.ll2004^ was designed mainly 
to rapidly detect, locate, and observe gamma-ray bursts 
(GRBs). After more than two years of operation, the Swift 
satellite has observed approximately 180 GRBs, for which 
more than 50 events have known redshifts reported in the 
Mission homepag^B- In this paper we present results from 
an homogeneous spectral analysis that we have carried out 
to the Swift long-duration GRBs with known z with the 
main aim of estimating the time-averaged spectral parame- 
ters of the prompt emission of these bursts, as well as related 
global quantities like the prompt energy fluence, S^, and the 
rest-frame isotropic-equivalent radiated energy, -Bt^iso. 

Unfortunately, the spectral coverage o f the Swift 7-ray 
detec tor, the Burst Alert Telescope (BAT. iBarthelmv et al.l 
I2OO5I ). is very narrow. The recommendation is to limit the 
spectral analysis to channels between 15 keV and 195 keV, 
though the Swift-BAT's team suggests a more conservative 
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upper limit of 150 keV because, above this energy, the cal- 
ibration is not yet sufficiently reliable. Thus, the BAT en- 
ergy range is not broad enough to aUow the spectra to be 
fitted unambiguously by the usual GRB broken power-law 
function of 4 parameters. Therefore, the inference of 7-ray 
spectral parameters and global quantities related to these 
parameters for the Swift GRBs is not an easy task. 

A main concern of this task is the determination of the 
uncertainties of the spectral parameters and their propaga- 
tion into the composite quantities. Due to the narrow energy 
range of Swift BAT we know a priori that the errors in the 
fitted spectral parameters will be large and correlated among 
them. Therefore, the appropriate handling of errors is cru- 
cial in order to make useful the spectral data generated by 
BAT. The effort is justifled by the great value that an homo- 
geneous sample has, where the spectral information for all 
the events is obtained with the same detector and analyzed 
with the same techniques and methods. This is the case of 
the growing-in-number sample of Swift GRBs. 

The characterization of the time-averaged photon 
spectra of GRBs with known redshift is important for 
several reasons. On one hand, this characterization offers 
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clues for understanding the radiation and particle accel- 
eration mechanisms at work during th e prompt phase of 
the bursts (see for r ecent reviews e.g ., JZhang fc MeszarosI 
I2OO4I : iPiran I l2005l : iMeszarod I2OO6I : Izhang 20o3). On 
the other hand, the spectral parameters and the global 



quantities inferred from them {S-,, 



etc.) are 



among the main properties that characterize GRBs. The 



other properties (e.g.. |Amati et al.l 12002 


: lAtteial l2003l: 


Ghirlanda. Ghisellini & Lazzatil 2004a: 


Yonetoku et al.l 


20041: Liang & Zhand 20051: Firmani et al.l 2006al) is cur- 



rently allowing to learn key aspects of the nature of GRBs 
(e.g.. [Thompson 112 006: Th ompson. Rees fc MeszarosI [2OO7I : 
iRamirez-Ruiz fcGranot ,20061 ). Furthermore, those corre- 
lations that are tight enough can be used to standardize 
the energetics of GRBs, making possible to apply GRBs 
for constructing the Hubble diagram up to unprecedentedly 
high redshifts (Ghirlanda et al. 2004b 2006; Dai et al . 2004; 
Firmani et al. 2005, 2006b, 2007; iLiang fc Zhangj l2005l : 
IXu. Dai fc Liand l2005l : ISchaeferl |2007^ . These tight corre- 
lations, for a given cosmological model, can be used also 
to estimate the (pseudo) redshifts of GRBs with not mea- 
sured redshifts, mai nly from the extensive CGfiO-BATSE 
databasjfl iLlovd-Ronning. Fryer fc Ramirez-Ruij 

l2002l : lAtteial liooa ). Thus, by combining the distribution 
in redshift of a large sample of GRBs with the observed 
flux distribution, inferences on the luminosit y function and 
formation rate of GRB s can be obtained (|Firmani et al.l 
l2004l : lGuetta et al.ll2005l) . 

The outline of this paper is as follows. In §2 we ex- 
plain the main steps of the spectral analysis carried out by 
us: selection of the sample (§§2.1), spectral deconvolution 
(fit) procedure and fitting models (§§2.2), and error analysis 
(§§2.3). The results of the spectral analysis are given in §3, 
along with the estimates of and -E^.iso as well as some 
correlations among the obtained quantities. The summary 
and conclusions of the paper are presented in §4. In the Ap- 
pendix the deconvolved spectra of the 47 GRBs studied here 
are plotted. 

We adopt the concordance ACDM cosmology with 
= 0.3, = 0.7, and Ho = 70 km s~^Mpc"\ 



2 SPECTRAL ANALYSIS 

The first step in our spectral analysis is the selection of the 
Swtft sub-sample of GRBs to be studied (§§2.1). Then we 
proceed to the spectral deconvolution with XSPEC and de- 
termination of the GRB spectral parameters for two- and 
three-parameter photon models (§§2.2). By using these pa- 
rameters, the energy fiuence and isotropic-equivalent energy 
are calculated. Finally, we construct the error confidence el- 
lipses for the spectral parameters and -B^^iso (§§2.3). 

The BAT is a large aperture 7— ray telescope. It consists 
of a (2.4xl.2)m coded aperture mask supported one meter 



In a recent paper, Q l l2007al ') showed that due both to a red- 
shift degeneracy in the so— called Amati correlation jAmati et al.l 
[2002) and to its high scatter, is not possible to infer reliable 
(pseudo)redshifts from this correlation. 



above the 5200 cm area detector plane. The detector con- 
tains 32,768 individual cadmium-zinc-telluride detector ele- 
ments, each one of (4x4x2)mm size. The effective efficiency 
of BAT starts at approximately 15 keV, attains a broad 
maximum at approximately 30-100 keV and then falls off, 
attaining at 195 keV the same level as at 15 keV. Due to the 
1 mm thickness of the lead tiles, the BAT coded mask starts 
to be transparent around 150 keV. We reduce the BAT's 
data using the analysis techniques provided by the BAT's 
instrument teairlf|. 



2.1 The sample 

We analyse the time-integrated spectra of all available 
Swtft GRBs with known redshift z. The list of events was 
taken from the daily updated compilation of J. GreineiQ 
for the period January 2005 to January 2007. The corre- 
sponding file spectra were loaded from the legacy NASA 
databasqj. The total number of events in this list was 
55. We exclude from the list the short-duration bursts 
(Tgo < 2s; GRB050509, GRB050724, and GRB061217) as 
well as those events with incomplete observational informa- 
tion (GRB050730, GRB060124, GRB060218, GRB060505 
and GRB060512), leaving us with a preliminary list of 47 
long GRBs with measured z. 

The list of 47 events is given in Table [T] The first four 
columns of Table [1] give respectively the burst name, the 
class (see below), the redshift z, and the time Tintthat we 
use to integrate the spectrum. In principle T90 could be used 
as an integration time. Nevertheless two problems arise in 
this case. The first one has to do with the systematic under- 
estimate of the fiuence, hence of E^^iao- The second problem 
may be more important; the elimination of the (early and 
late) wings in the light curve can influence the estimates of 
the spectral parameters. For example the elimination of the 
late wing may harden artificially the spectrum. For these 
reasons we have fixed the integration time by a visual in- 
spection of the light curve as the period where the signal 
is clearly identifiable in the background noise. A modest 
overestimate of Tint does not influence on the values of 
the integrated parameters, though increases the noise. We 
would like to remark that Tint is not a good estimate of 
the burst duration because its uncertainty is related to the 
identification of faint (some time extended) wings. 



2.2 Spectral fits 

It is well known that the GRB photon spectra are in general 
well described, in the energy range of ~ 10 keV to a few 
MeV, by the so-called Band function (|Band et al.lll993l ). 
which is a two smoothly connected power laws: 



N{E) = iVo ( 



E 



100 keV 



iVo 



So(q-/3) 



100 keV 



{a-/3) 



E 

E < Eo{a 
E 



(1) 



^^P^^-^Kiootev) 



^ http:/ /heasarc. gsfc.nasa.gov/docs/swift/analysis/threads 
* http://www.mpe.mpg.de/ jcg/grbgen.html 
^ ftp:/ /legacy.gsfc. nasa.gov/swift/data/obs 



Spectral analysis of Swift GRBs with known z 3 



GUB Class z Tint LgiVo[sd] a[sd\ ei m 62 na x?cc 

see 50 keV 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



050126 


A 


1.29 


050223 


A 


0.5915 


050315 




1.949 


050318 


A 


1.44 


05031 q 




3.24 


050401 


B 


2.9 


05041 6 




6535 


050505 


A 


4.27 


050525 


A 


0.606 


050603 


B 


2.821 


050803 


A 


0.422 


050814 


A 


5.3 


050820 




2.61 


050824 




0.83 


050904 




6.29 


050908 


A 


3.344 


050q22C 


B 


2.198 


051016B 


D 


0.9364 


051109 


A 


2.346 


051111 


B 


1.55 


051221 




0.547 


060115 


A 


3.53 


060206 


A 


4.048 


060210 


A 


3.91 


060223 


A 


4.41 


060418 


c 


1.489 


060502 A 


B 


1.51 


060510B 


A 


4.9 


060522 


A 


5.11 


060526 




3.21 


060604 




2.68 


060605 


A 


3.711 


uuuuu / 


A 


o.uoz 


060614 


C 


0.125 


060707 


A 


3.425 


060714 


A 


2.71 


060729 


C 


0.54 


060904B 


A 


0.703 


060906 


A 


3.686 


060908 


A 


2.43 


060926 


D 


3.2 


060927 


A 


5.6 


061007 


C 


1.261 


061110B 


C 


3.44 


061121 


C 


1.314 


061222B 


A 


3.355 


070110 


C 


2.352 



40.0 


-2.656(0.026] 


-1 


30.0 


-2.652(0.028] 


-1 


129.0 


-2.567(0.013] 


-2 


34.5 


-2.375(0.014] 


-1 


190.0 


-3.100(0.034] 


-2 


43.0 


-1.667(0.012] 


-1 


5.0 


-2.225(0.055] 


-3 


70.0 


-2.418(0.018] 


-1 


15.0 


-0 932f0 0031 


_1 


16.0 


-1.347(0.010] 


-1 


165.0 


-2.868(0.022] 


-1 


95.0 


-2.750(0.030] 


-1 


32.0 


-2.461(0.016] 


-1 


38.0 


-3.102(0.081] 


-2 


207.0 


-2.603(0.009] 


-1 


24.0 


-2.685(0.031] 


-1 


10.0 


-1.761(0.009] 


-1 


6.0 


-2.509(0.048] 


-2 


43.0 


-2.052(0.033] 


-1 


45.0 


-2.080(0.008] 


-1 


3.1 


-1.422(0.008] 


-1 


128.0 


-2.851(0.018] 


-1 


14.0 


-2.184(0.011] 


-1 


278.0 


-2.579(0.015] 


-1 


14.8 


-2.320(0.020] 


-1 


80.0 


-1.993(0.006] 


-1 


34.0 


-2.157(0.011] 


-1 


289.0 


-2.839(0.011] 


-1 


80.0 


-2.829(0.024] 


-1 


320 


-3.387(0.035] 


-1 


67.0 


-3.356(0.085] 


-2 


37.0 


-2.810(0.026] 


-1 


127.0 


-2.686(0.011] 


-1 


154.0 


-1.877(0.005] 


-1 


81.0 


-2.652(0.022] 


-1 


135.0 


-2.614(0.016] 


-1 


134.0 


-2.658(0.021] 


-1 


187.0 


-3.059(0.022] 


-1 


63.0 


-2.386(0.020] 


-2 


30.0 


-2.040(0.009] 


-1 


13.0 


-2.724(0.044] 


-2 


32.0 


-2.388(0.013] 


-1 


85.0 


-1.334(0.004] 


-1 


147.0 


-3.055(0.021] 


-1 


125.0 


-1.936(0.004] 


-1 


57.0 


-2.367(0.023] 


-1 


103.0 


-2.806(0.016] 


-1 



317 


Q 


lOll 





026 


_l 


00 





01 


837 


Q 


1061 





025 


-0 


99 





13 


101 


Q 


0551 





012 


-0 


99 





11 


968 


Q 


0571 





012 


-0 


99 





13 


345 


Q 


1181 





026 


-0 


98 





17 


472 


Q 


0461 





Oil 


_1 


00 


-0 


04 


119 


Q 


1551 





026 


-0 


95 





30 


415 


Q 


0711 





018 


_1 


00 


-0 


03 


769 


Q 


0151 





003 


_l 


00 


-0 


09 


159 


Q 


0431 





010 


_1 


00 


-0 


05 


526 


Q 


0851 





021 


_l 


00 





04 


798 


Q 


1081 





027 


-0 


99 





12 


689 


Q 


0671 





015 


_l 


00 





09 


743 


Q 


2191 





036 


-0 


95 





30 


205 


Q 


0411 





009 


_1 


00 


-0 


04 


947 


Q 


1071 





025 


-0 


99 





17 


371 


Q 


0381 





009 


_l 


00 


-0 


02 


331 


Q 


1451 





030 


-0 


97 





24 


626 


Q 


1241 





032 


_1 


00 





05 


305 





034] 





008 


_l 


00 


-0 


04 


397 


Q 


037] 





008 


_l 


00 


-0 


03 


771 





071] 





016 


-0 


99 





11 


693 


Q 


047] 





Oil 


_l 


00 





06 


544 


Q 


059] 





015 


_l 


00 


-0 


01 


775 


Q 


078] 





019 


-0 


99 







634 


Q 


030] 





006 


_1 


00 


-0 


01 


431 


Q 


045] 





Oil 


_1 


00 


-0 


01 


778 





048] 





010 


_l 


00 





07 


578 


Q 


093] 





023 


_l 


00 





08 


924 


Q 


140] 





028 


-0 


99 





17 


340 


Q 


283] 





060 


-0 


98 





21 


474 


Q 


105] 





025 


_1 


00 





07 


4oD 


U 


U4DJ 


u 


UiU 


1 
i 


UU, 


u 


no 


936 





024] 





004 


-1 


00, 





01 


667 





082] 





021 


-1 


00, 





06 


935 





063] 





013 


-0 


99, 





12 


870 





083] 





019 


-0 


99, 





12 


628 





090] 





021 


-1 


00, 





10 


006 





075] 





017 


-0 


99, 





14 


339 





039] 





009 


-1 


00, 


-0 


01 


460 





146] 





026 


-0 


97 





24 


654 





052] 





012 


-1 


00 





07 


002 





016] 





003 


-0 


99, 


-0 


13 


069 





091] 





021 


-1 


00, 


-0 


05 


410 





018] 





004 


-1 


00, 


-0 


08 


962 





084] 





020 


-0 


99, 





14 


581 





702 





016 


-1 


00, 





06 






106 


-0 


01 


-1 


00) 


1 


399 





109 


-0 


13 


-0 


99) 





902 





056 


-0 


11 


-0 


99) 





889 





060 


-0 


13 


-0 


99) 





936 





119 


-0 


17 


-0 


98) 





687 





047 





04 


-1 


00) 


1 


060 





182 


-0 


30 


-0 


95) 





971 





073 





03 


-1 


00) 





789 





015 





09 


-1 


00) 


2 


934 





043 





05 


-1 


00) 





913 





088 


-0 


04 


-1 


00) 





882 





117 


-0 


12 


-0 


99) 


1 


036 





067 


-0 


09 


-1 


00) 





988 





241 


-0 


30 


-0 


95) 


1.0157 





041 





04 


-1 


00) 


1 


070 





114 


-0 


17 


-0 


99) 





969 





038 





02 


-1 


00) 





900 





158 


-0 


24 


-0 


97) 


1 


344 





129 


-0 


05 


-1 


00) 





819 





035 





04 


-1 


00) 





764 





037 





03 


-1 


00) 





996 





073 


-0 


11 


-0 


99) 





875 





049 


-0 


06 


-1 


00) 


1 


061 





060 





01 


-1 


00) 





928 





082 


-0 


11 


-0 


99) 


1 


009 





030 





01 


-1 


00) 





815 





046 





01 


-1 


00) 


1 


031 





049 


-0 


07 


-1 


00) 





643 





098 


-0 


08 


-1 


00) 


1 


188 





132 


-0 


17 


-0 


99) 





967 





290 


-0 


21 


-0 


98) 





804 





109 


-0 


07 


-1 


00) 





912 





047 


-0 


02 


-1 


00) 





947 





024 





01 


-1 


00) 


1 


072 





087 


-0 


06 


-1 


00) 


1 


116 





064 


-0 


12 





99) 


1 


161 





084 


-0 


12 


-0 


99) 





863 





090 


-0 


10 


-1 


00) 





855 





078 


-0 


14 


-0 


99) 


1 


114 





040 





01 


-1 


00) 





864 





151 


-0 


24 


-0 


97) 





972 





055 


-0 


07 


-1 


00) 


1 


158 





016 





13 


-0 


99) 





523 





091 





05 


-1 


00) 


1 


175 





018 





08 


-1 


00) 





450 





089 


-0 


14 


-0 


99) 


1 


280 





070 


-0 


06 


-1 


00) 


1 


007 



Table 1. The sample of Swift long GRBs with known z analysed here (47 events). The z and duration (Tint) of the burst as well as the 
results from the PL photon model fit are reported. Second column refers to the class group assigned to the burst according to our spectral 
analysis (see text). The best fit values of the parameters and their one— dimensional standard deviations (sd) are given in columns 5 and 
6; 61^2 and ni 2 in columns 7 and 8 are the semi-axis lengths and principal axes of the error ellipse, respectively. The reduced 
reported in the last column (dof=58). 



E> Eo{a~i3) 

where a and /3 are the photon indices of the low and high 
energy power laws, respectively, and Eq is the e— folding 
(break) energy related to the peak energy in the vf,^ [or 
E'^f{E)] spectrum by Epy, = Eox{2 + a). Note that £pk is 
well defined for a ^ — 2 and /3 < — 2. The normalization A'o 
is in photons s"'^ cm~^ keV~^. 



The spectral range of BAT is narrower than that one 
of BATSE and other previous detectors. In the most recent 
update of the BAT -team homepag^fl it is written that due 
to varying threshold levels in individual detectors, channels 
below 15 keV should not be used for spectral analysis. Like- 

^ http: / /swift. gsfc.nasa.gov/docs/swift/analysis/bat-digest. html 
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wise, channels above 150 keV are unreliable due to a lack of 
calibration data at those energies. Then, we are limited to 
the narrow energy band of 15-150 keV for the spectral anal- 
ysis of Swift bursts. This implies that model fitting with the 
4-parameter Band function will not be reliable for most of 
these burst. Therefore, we proceed to fit the observed pho- 
ton spectra to both the power law (PL) and cutoff power 
law (CPL) models. The former has only two parameters, the 
normalization, A^o, and the photon index, a, and is given by: 



(2) 



The CPL model implies one more parameter, an energy 
e— folding spectral break, Eq, related to _Epk in the same 
way as in the Band model (see above). The CPL model is 
defined as: 



NiE) = iVo 



E 



. 50 keV 



a 


■ E ■ 


1 exp 





(3) 



It is easy to verify that the CPL model (named also the 
Compton -COMP- model) is the Band model with /3 
—oo. 

As we will see later, the complex nature of the errors 
due to the small energy range of the spectra makes it more 
convenient to use the logarithms of A'o and Eq (or -Epk) 
rather than the linear quantities. This choice reduces the 
asymmetry of the errors and of the confidence levels (CL) 
contour shapes. Therefore, we have introduced into XSPEC 
new models corresponding to Eqs. ([2]) and Q (and also Eq. 
(1); see §§3.3) to carry out the fits with LogA^o and LogSpk. 

The spectral analysis was carried out by using the hea- 
soft6.1.2 public softwar^ll- The first step in our analysis is to 
define the background-subtracted light curve from the cor- 
responding data file. The photon counts at each channel are 
taken in the time interval Tint where the main light curve 
is clearly above the background noise. The spectral file was 
corrected by position and systematic errors. Then, the pho- 
ton counts at each channel is convolved with a response ma- 
trix, build up for a given event with the calibration matrij|f|, 
in order to obtain the time-integrated energy spectrum. As 
mentioned above, the BAT detector is well calibrated only 
in the 15-150 keV energy range. Therefore, we used the 60 
BAT channels in this range, with the default binning. The 
spectral fit to the different photon models mentioned above 
was performed with the heasoft package XPSEC (version 
11.3.1). 

For the GRBs for which the (CPL) _Epk may be esti- 
mated, we calculate the fluence, 5*^., corresponding to the 
observed spectral band (15-150 keV) only with the aim of 
completeness. In addition the rest isotropic-equivalent en- 
ergy, Ery^iso, is calculated from the observer bolometric fiu- 
ence S'y ^ corresponding to the spectral range 1-10000 keV 
at rest: 



E^ 



4TvS-ygdLizf 
(1+^) 



(4) 



where dL{z) is the cosmology-dependent luminosity dis- 
tance. 



http: / /heasarc. gsfc.nasa.gov/docs/software/lheasoft/download. html 
* We use the January 2007 matrix calibration provided by the 
Swift-BAT team. 



2.3 Errors from the fit and error propagation 

A main concern of the spectral deconvolution of BAT spec- 
tra is the expected large uncertainty and correlations among 
the fitted spectral parameters. Consequently, a careful inter- 
pretation of the errors is crucial to make usable and reliable 
the spectral analysis of Swift GRBs. 

First, as mentioned above, we have carried out the fits 
using LogA'^o (PL model), and LogA'o and LogiSpk (CPL 
model), instead of the corresponding linear values. The er- 
rors in the linear parameters are very asymmetric and corre- 
lated among them in a very complex way; the corresponding 
regions at a given CL, as plotted with the XSPEC com- 
mand steppar in the plane of two parameters, have shapes 
strongly deviating from ellipses. Instead, the errors and the 
CL regions when the fits are carried out with logarithmic 
quantities, tend to be more symmetric and ellipse-shaped, 
respectively. This facilitates the handling of the errors to 
calculate composite quantities and correlations among the 
GRB parameters. 

The fit command of XSPEC, generates ellipsoidal con- 
fidence intervals that approximate the complex regions of 
joint variation of the parameters. These ellipsoids are char- 
acterized by the semi-axis lengths, (they are the square 
root of the variances), where i is for each axis, and the cor- 
responding principal axes (unitary vectors rii). We use these 
data to calculate the covariance matrix, dj, and the ellip- 
tical boundary of the desired CL region in the subspace of 
interest (see e.g.. [Press et al. 1999). The diagonal elements 
of the matrix, Cu, are the variances of each model param- 
eter (two for the PL model and three for the CPL model), 
while the off-diagonal elements (covariances) show the cor- 
relations among the parameters. 

The error propagation from the spectral quantities 
(LogA'^o, a, Log_Epk) to the composite quantities LogiSpk 
(being E'^y. the peak energy at rest _Epk = Ep^x [l-|-z)]) 
and Log_E-y4so is performed in two steps. First, we have ex- 
tracted sets of spectral quantities by a Monte Carlo method 
according to their elliptical CL. Then, for each set we have 
calculated the corresponding composite quantities as well 
as their covariance matrix averaging on the overall extrac- 
tions. From this last covariance matrix the elliptical CL of 
the composite quantities is obtained. 



3 RESULTS 

In Table [T] we present the results from the spectral fit to the 
PL model for the whole sample of 47 long GRBs. Columns 
1 to 4 give the burst name, the class group of the burst, z, 
and Tint. Columns 5 and 6 give the best fit values of the 
normalization, LogA'^o, and the photon index, a, with the 
corresponding one-dimensional symmetric standard devia- 
tions in brackets, as obtained with the XSPEC error com- 
mand. In columns 7 and 8 we report the square root of the 
variances, e^, together with the unitary vectors ni (principal 
axes), which provide respectively the semi-axis lenghts and 
directions of each axis (i = 1,2) of the joint error ellipse. 
The reduced is reported in column 9. 

As it may be appreciated from Table [TJ the best-fit a is 
smaller than -2 for 7 bursts out of the 47 (GRB 060906 is 
marginal and it will be recovered by the CPL fitting) . Most 
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Figure 1. Error CL contours projected in the plane of the pair of variables Logi?pk-(-a) for each one of the bursts in samples A and 
B (CPL model). The irregular contours (black lines in the electronic version) were calculated with the XSPEC command steppar, while 
the elliptical contours (red lines) were constructed from the variances and principal axes given by the XSPEC command fit. Inner (thick 
line) and outer (thin line) contours in each case are for Ax^ = 1 and 2.3, respectively (see text). The dotted line indicates the level of 
150 keV. 
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Table 2. Spectral fit results using the CPL photon model for the class A and B GRBs. The best fit values of the parameters and 
their one-dimensional standard deviations (sd) are given in columns 2-4; ei_2,3 and ni_2,3 in columns 5-7 are the semi-axis lengths and 
principal axes of the error ellipsoid, respectively. The reduced is reported in the last column (dof=57). 



likely, the vf^ spectra of these GRBs peak at energies below 
15 keV. They could be related to the X-ray rich events and 
X-ray Flashes discovered by BeppoSAX and HETE-2 (Kip- 
pen et al. 2001; Sakamoto et al. 2005). 

We have also fitted the 47 GRB spectra to the CPL 
model and found that for 18 bursts the fit produces too 
large uncertainties and/or unphysical solutions. For these 
18 cases the PL model is definitively more reliable than the 
CPL one. According to the results from our spectral fits 
(judged by the joint 68% CL contour errors), we classify the 
47 Swift GRBs analysed here in four classes: 

Class A - The CPL-model fit is acceptable and -Epk is 
within the 15-150 keV range. 

Class B - The CPL-model fit is acceptable, but iJpk > 
150 keV. 

Class C - the CPL-model fit gives unreliable results, 
while the PL-model fit is acceptable and gives a> —2, sug- 
gesting that i?pk is at energies > 150 keV. 

Class D - The CPL-model fit gives unreliable results, 
while the PL-model fit is acceptable and gives a< —2, sug- 
gesting that Spk is smaller than 15 keV. 

The class of each GRB is indicated in the second col- 
umn of Table [1] This classification can be easily understood 
after a visual inspection of the plots shown in Appendix and 
figures [T]and[l In Figs. EH |H and |M] of the Appendix we 
present the deconvolved time-integrated vf^ [or E^f{E)] 



observed spectra of the 47 Swift GRBs analysed here. The 
CPL model was used for the spectral deconvolution. The 
error bars show the observations while the continuous lines 
show the fitted curves. The spectra of all the events that 
we classify in groups A and B show some curvature, evi- 
dence of a peak in the f{E) distribution. Nevertheless, 
in most of the cases the formal f — test analysis says that 
the spectral fits from the PL model to the CPL one (which 
introduces one more parameter, i5pk) do not improve sig- 
nificantly. Therefore, one expects large uncertainties in the 
determination of i?pk for most of the events analyzed here. 
This is why it is important to carry out a carefully statistical 
analysis of the errors and to adequately propagate them for 
calculating composite quantities and inferring correlations. 

Figures [T] and [5] show different contours of the CLs pro- 
jected in the plane of the pair of variables LogEpi^ vs. (—a) 
for each burst in groups A and B (CPL model). The irregu- 
lar CL contours (black lines in the electronic version) were 
calculated with the XSPEC command steppar, while the el- 
liptical CLs (red lines) were constructed from the joint vari- 
ances and principal axes given by the XSPEC command fit 
(see §§2.3; columns 5, 6, and 7 in Table [2] below). Thick lines 
show the CLs for Ax^ = 1, while thin lines show the CLs for 
Ax^ = 2.3. For Ax^ = 1, the projections of the CL regions 
in each axis enclose « 68% of the one-dimensional probabil- 
ity for the given parameter (one standard deviation), while 
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Figure 2. Continues Fig. [T] 



for Ax^ = 2.3, the plotted region encompasses « 68% of 
the probabiUty for the joint va riation of the two parameters 
(one standard deviation) (e.g.. |Press et al.]|l999l ). As can be 
appreciated, the eUipses are a reasonable description of the 
CL regions in most cases, with the great advantage that they 
are handled analytically. 

We see from Figs. [1] and [2] that for many of the GRBs, 
the LogiJpk and a spectral parameters are correlated (in- 



clined, thin ellipses). A similar situation arises for the other 
combination of parameters. Therefore, for this kind of data, 
a simple statistical quantification of the errors, such as the 
one-dimensional standard deviation, is not enough. This is 
why we have worked out the full covariance matrix calcu- 
lated from the error ellipse. We use it to propagate errors in 
the calculation of E'^^^ and -E^^iso (§§2.3). This is a crucial 
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Table 3. Calculated fluence and rest— frame isotropic and peak energies for our GRB groups A and B. The corresponding one-dimensional 
standard deviations are given between square brackets. The last 2 columns give the estimated semi— axis lengths (ei 2) and principal 
axes (ni^2) of the error ellipse associated to the pair of quantities LogE^j^ and Logi?.y_iso. Note that the errors in most of the cases are 
correlated significantly (see also the plotted ellipses in Figs. 4 and 5). 



step when using the data, for instance, for establishing cor- 
relations among the spectral or composite quantities (§§3.3). 

The CPL-model fittings for GRB groups A and B are 
presented in Table Column 1 gives tlie burst names. 
Columns 2, 3, and 4 report tlie best fit values of LogA'o, 
a, and LogiSpk with their corresponding one-dimensional 
standard deviations, as calculated with the XSPEC error 
command. In columns 5, 6, and 7 we report the lengths 
ei of each semi-axis together with the unitary vectors rii 
(principal axes), which provide the directions of each axis 
[i = 1, 2, 3). The reduced is reported in column 8 (dof — 
57). 

Finally, in Table [3] we report the calculated values for 
S^, Logi?^,iso, and LogSp^ , as well as the corresponding 
propagated errors for the samples A and B, using the CPL- 
model fit results (see §2.3 for the used error propagation 
procedure). Column 1 gives the burst name. Columns 2 and 
3 give respectively the values of (in the observed range 
15-150 keV) and of S-y.iao (extrapolated to the rest-frame 
energy range 1-10000 keV) corresponding to the best fit 
spectral parameters. Column 4 gives the value of the rest- 
frame peak energy, E'^^. The one-dimensional standard de- 
viationss of Sj, iJ-y, iao, and E'^y. are given within square 
brackets. However, due to the high correlation among the 
errors the simple standard deviation is not enough to char- 



acterize statistically the errors. Thus, in columns 5 and 6 
we report the error ellipse elements (semi-axis lenghts, e^, 
and principal axes ni, i — 1,2) associated to the pair of 
quantities LogEp^. and LogiJ.^,iao. It is important to remark 
that the correlation between the uncertainties of hogE'^^ 
and Logi?7,iso is not entirely due to the correlation among 
LogA^o, ct, and LogSpk. For instance an increase in LogiSpk 
alone is able to produce an increase in both LogE'^j^ and 
Log£-y4so. This effect will turn out to be helpful in the next 
section, where we will study how hogE'py. and Log_E.y,iso do 
correlate in the entire GRB sample. 

The median and dispersion (quartiles) of the best- 
fit CPL parameter a for the 24 + 5 classes A+B 
GRBs are — 1.12]|^q j^^, whicli is i n agreement w ith the re- 
suhs from BAT SE bright-G RB ('Kancko et al.' '2006') and 
HETE-2 GRB ()Barraud et a l. 2003; Sakamoto ct al. 200i) 
time-integrated spectra. Our a distribution is in particular 
similar to the one from HETE-2, with only ~ 25% of the 
bursts having a < —1.2. The median and quartiles of the 
LogBpk [keV] parameter are 1.91+0 Unfortunately, due 
to the low energy limit of the Swift-'BAT, the sample anal- 
ysed here covers only the ~ 15% fraction of bursts in the 
(CPL-i nodel) _E;pk distribu tion inferred for BATSE bright 
bursts ()Kaneko et al.ll2006l ). 

The distribution of E'^y. (at rest) is much broader than 
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Figure 3. Top panel: the observed i?p]j as derived by tlie CPL fit 
as a function of the spectral index acPL of the same fit. Bottom 
panel: iJpj; (as derived by the CPL fit) as a function of the spectral 
index apL of the single power law fit. The correlation of i?p]j with 
apL is induced by i?p]j entering in the Swift energy band (see 
text), and cannot be used to infer Ep]^ when knowing only apL- 

the one of i?pk (the median and quartiles of Log^Epj^ [keV] are 
2.54lg 21). This is expected due to the broad distribution in 
redshfits of the Swift bursts analysed here. For our 24 + 5 
classes A+B GRBs, the z distribution has a broad maximum 
at 2; ~ 3 — 4; the median and quartiles are S.OSlj'g^. The 
median and quartiles of LogB^^iso [10^" erg] for our 24 + 5 
classes A+B GRBs are 2.78lo;44. For the 8 long-duration 
Swift GRBs in common with the compilation presented in 
Amati (2006b), the values of E'^y. reported here and in that 
paper are similar, but those of E^.iso are systematically 
larger in Amati (2006b). The latter difference is probably 
due to the fact that Amati (2006b) did not use for most of 
those GRBs the spectral information from the Swift satel- 
lite but from other space missions (Konus, HETE-2), which 
allowed to find a fitting with the Band model; as will be 
discussed in §§3.3, this implies typically a larger E^^iso with 
respect to the case a CPL model is used. 



Our analysis with the Swift data (see Tab. 3) for the time 
integrated spectrum and a CPL model yielded a = —1.23 ± 
0.19 and Spk = 186^7^ keV. Within the (relatively large) 
errors, the results are consistent. 

The Konus-Wind data of GRB 050603 have been fitted 
by Golenetskii et al. (2005b) with a Band model in the [20- 
3000] keV energy range, with a = -0.79 + 0.06, /? = -2.15 + 
0.09 and Bpk = 349 + 28 keV. Our analysis with a CPL 
model gives a = -1.23 + 0.19 and £pk = 323^^64 keV. The 
derived low energy spectral index is somewhat softer, while 
the value of _Epk is consistent. 

GRB 050922C was detected by HETE-2 with the FRE- 
GATE instruments in the 7-30, 7-80, and 30-400 keV bands. 
It was fitted with a CPL law model by Crew et al. (2005) 
with a = -0.83lo;26 and E'pk = 130.5t26:8 keV, to be com- 
pared with our values of a = -1.16+0.15 and £pk = 224t|f 
keV. Again, within the errors, there is consistency. 

GRB 051109, detected by Konus-Wind and fitted with 
a CPL model by Golenetskii et al. (2005c) in the [20-500] 
keV energy range, showed a = — 1.25lQ5g ^^'^ ~ 
16ll5g* keV. This should be compared with our a = —1.04+ 
0.55 and Spk = 83+^7 keV. Even in this case there is con- 
sistency, within the relatively large errors. 



3.2 Peak energy vs spectral index 

In Fig. Owe show -Epk, derived with the CPL fit, as a func- 
tion of the spectral index acPL of the CPL fit (top panel) 
and as a function of apL, the index derived when fitting 
with a single power law. As can be seen, although there is 
no relation between iSpk and qcpl, a correlation appears 
when fitting with a single power law. This correlation, how- 
ever, has no physical meaning, since it is the result of the 
attempt of the single power law to account for the data 
point above the peak, with smaller flux (in ufi, ). As -Epk 
decreases, a larger fraction of data lies above the peak, in- 
ducing the single power law to steepen. This also means 
that this effect vanishes when _Epk is outside the Swift en- 
ergy range, therefore if we do not know -Epk, we cannot use 
this correlation to infer it. Furthermore, note the difference 
in the derived values of the spectral indices. As expected for 
these bursts for which we could derive -Epk, the index apL 
is systematically softer than acPL . Given the obvious result 
of this subsection, we have considered unnecessary to show 
the upper panel of Fig. |3]with the 68% CL error ellipses. A 
complete treatment of the errors in this sense confirms the 
previous result. 



3.1 Bursts detected also by other instruments 

Among the bursts listed in Tables 2 and 3 there are four 
GRBs which have been detected by other instruments be- 
sides Swtft. Three of them have been detected by Konus- 
Wind (GRB 050401, 050603 and 051109A) and one by 
HETE-2 (GRB 050922C). 

For GRB 050401 the Konus-Wind results reported by 
Golenetskii et al. (2005a) concern the spectral parameters 
for the two peaks displayed by its light curve, fitted with 
a Band model in the energy range [20-2000] keV. The first 
peak had a = -1.15 + 0.16, /? = -2.65 + 0.31 and E^pk = 
132 + 16 keV, while the second peak is described by a = 
-0.83 + 0.21, P = -2.37 + 0.14 and Epk = 119 + 26 keV. 



3.3 Correlation between E'^y. and S^^so 

In Fig. |4l -Epi; is plotted vs -E-yjso for the class A (red-line 
ellipses) and class B (green-line ellipses) Swift GRBs. The 
ellipses correspond to the joint 68% CL error regions calcu- 
lated as explained in §§2.3. The dotted straight lines encom- 
pass the 3cr scatter of the updated Amati correlation pre- 
sented in iGhirlanda I (|2007| ) for 49 GRBs; the central thin 
solid straight line is the best fit for that sample. A visual 
inspection shows that the Swift data analyzed in this paper 
are indeed correlated in the -Epi^--E^,iso diagram. Taking 
into account the uncertainties, no GRB may be classified 
as outlier. Eight GRBs show their central point above the 
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Figure 4. Correlation between E'^^ and -E-,jso (rest frame) as obtained from spectral fits using the CPL model for the class A (red 
ellipses) and B (green ellipses) bursts. The ellipses correspond to the jo int 68% CL error regions of each data point. The dotted lines 
encompass the Scr scatter of the updated Amati correlation presented in lGhirlanda I jSooJ) for 49 GRBs; the thin solid line is the best 
fit for that sample. The thick solid line is the best fit that we find for the Swift data presented here taking into account the correlated 
errors, and the dot-dashed line is for the la uncertainty, computed in the barycentre of the data points. 



3(7 strip, but taking into account the 68% CL contours as 
well as the fact that the spectral fitting was performed us- 
ing a CPL model instead of the Band one (see below), make 
such GRBs compatible with the 2>o strip. Notice that the 
5 (marginal) class B events are not outliers. Therefore, we 
may include them for a fit along with the 24 class A events. 

The thick solid line in Fig. |4] is the best fit that we 
find for the 29 class A and B GRBs. The linear fit has been 
performed by an iterative minimum method, where the 
coordinate origin of the data is shifted to the barycentre, and 
the uncertainty component of each datapoint is calculated 
from the corresponding error ellipse in the orthogonal direc- 
tion of the straight fitted line. The dot-dashed line curves 
describe the standard error of the mean values of the fit 
parameters computed in the barycentre of the data points. 
The best fit in logarithmic quantities is: 

= (2.25±0.01)+(0.54±0.02)Log {jg^^ (5) 



We point out the weight of GRB050525 in determining the 
best fit due to its very-reduced spectral parameter uncer- 
tainties. We have also carried out the linear fit for only the 
24 class A GRBs: 

(§7) = (2.23±0.01)+(0.54±0.02)Log (y^fei^) (6) 

As expected, the fit is very similar to the one that includes 
the 5 class B GRBs. 

Thus, a correlation between E'^]^ and £7,130 is confirmed 
for 29 long GRBs observed by the Swift satellite and ho- 
mogeneously analyzed by us. The values of the power-law 
index {a ~ 0.54) and normalization {K ~ 100 keV for 
i?7,iso given in units of lO^^erg, E'^^ — KE-y^iBo°') corre- 
sponding to the best fit, lie in the ranges of values found 
with different previous samples (see for a review of previ- 
ous studies in Amati 2006a). It is encouraging that with 
the Swift data, as analysed here, one obtains an E'p^^-E^,iao 
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correlation fully consistent with the one obtained for GRBs 
observed wit h other satellite d etectors (the so-called Am- 
ati relation; lAmati et al ] I2OOI see also Amati 2006a and 
Ghirlanda 2007 for the most recent updates), which had a 
much broader energy (spectral) range than BAT. The er- 
rors in the fit parameters reported in eqs. ^ and ([6]) are at 
the 68% CL. These errors are somewhat larger than those 
reported in Ghirlanda (2007) and Amati (2006a; note that 
therein the errors are at the 90% CL). 

In general, the scatter around the correlation obtained 
here is within the range of different previous determinations 
of the Amati correlation. The value of the reduced for 
the fit with the 29 class A+B GRBs is xl = 379/27 = 14 
close to the values xl = 493/41 = 12 and xl = 530/47 = 
11.3 obtained for larger samples (43 and 49 GRBs) by 
Ghirlanda et al. (2006) and Ghirlanda (2007), respectively. 
We have also estimated the v ariance of our fit, (see e.g., 
iBevington &: Robinson I Il992l . Chapter 11), by taking into 
account the joint variance of the data, i.e. the elliptical CL 
contours. For the fit using the 24+5 GRBs, we find s = 0.13 
(dof=27). Note that for calculating the deviations, we use 
orthogonal distances to the best-fitting line rather than dis- 
tances projected in the LogE'^^ (Y) axis. An estimate of the 
vertical la spread of the data around the best-fitting model 
would be j^i = 0.15, in agreement with previous esti- 
mates (see e.g.. lAmatill2006al ). 

The virtue of the Swift sample analysed here is its homo- 
geneity, in the sense of observations, reduction, and analysis. 
Previous samples used to infer the E'py-E^,iso correlation 
were constructed from different detectors, with an informa- 
tion collected from different authors, who used different data 
processing and assumptions. This heterogeneity in the data 
should introduce a spurious scatter. In spite that the data 
used here have large observational errors, the dispersion of 
our fit lies within the range of fits reported previously. This 
is namely due to the homogeneity of the sample and be- 
cause, having taken into account the whole error ellipses 
associated to the data points, it results that the major-axis 
orientation of most of these ellipses is close to the direction 
of the correlation line. 

Examining in more detail Fig. |31 one sees that the nor- 
malization of the correlation found here [eq. ([SJ] is slightly 
higher (by ~ 0.1 dex at the sample median value of iJ^.iso) 
than the previously established one. This is indeed expected. 
Due to the narrow energy range of BAT, we had to force 
the GRB spectra to be fitted by the GPL model, while most 
spectra could have been better described by the Band model, 
if higher energies had been detected. It is known that when 
fitting the CPL model to a spectrum whose shape is ac- 
tually described by the Band model, then (i) Epy^ is larger 
than it would otherwise be, and (ii) the low-energy power- 
law index a is m ore negative than it would otherwise be 
IIBand et al. l ll993l : lKaneko et al.ll200^ . On the other hand, 
(iii) in the CPL model the vfi, spectrum is cut-off exponen- 
tially after the break energy Eq, while in the Band model 
the spectrum decreases as a power law. Item (i) implies di- 
rectly that i?pk is likely overestimated when the CPL model 
is used in the spectral analysis, and items (ii) and (iii) im- 
ply that Er^^iao, which is calculated in the rest-frame energy 
range of 1-10000 keV by extrapolating the fitted spectrum, 
is likely underestimated. 



In order to explore in a statistical sense the effect in 
the E'p^^~E^,iso correlation of using the CPL model instead 
of the Band one for the spectral fitting, we carried out the 
following experiment. The time-averaged spectra of the 29 
long GRBs that had a reasonable CPL-model fit (groups A 
and B) were fitted with XSPEC to a Band model, but with 
the high-energy power law index /3 frozen to the typical 
value of —2.3. As in the case of CPL, we have modified the 
Band model in order to perform the fit for the logarithms 
of the normalization and Spk. Then, we followed the same 
procedures described in §2 for handling the errors, taking 
into account the correlation among them. Fig. [S] shows the 
final data points and their 68% CL ellipses in the -Epk-i?7,iao 
diagram. As expected, the normalization of the E'^^-Ery^iao 
correlation decreases (with respect to Fig. U)), now being 
al most equal t o the one of the updated Amati relation given 
in I Ghirla nda (2007) for the range of iS^.iso values studied 
here. The best fit in logarithmic quantities is: 

(§7) = (2-23±0.01)+(0.53±0.04)Log (y^ll^) (7) 



4 SUMMARY AND CONCLUSIONS 

We have analysed homogeneously the time-integrated spec- 
tra of a relatively large sample of long-duration GRBs with 
known redshift obtained by a same detector, the Swift-BAT 
instrument. We have carried out spectral fits to 47 GRBs 
by using the two-parameter PL and three-parameter CPL 
photon models. Due to the narrow spectral range of BAT, 
15-150 keV, only two- or three-parameter models can give 
reliable fits. However, even in these cases, large uncertain- 
ties are expected in the fitted parameters, making the study 
of the (highly correlated) errors important. We have per- 
formed Monte Carlo simulations in order to propagate the 
errors and calculate composite quantities such as E'^y. and 
i?7,iso. The main results and conclusions are as follow: 

• For 29 GRBs (classes A-l-B) out of the 47 bursts anal- 
ysed here, the CPL model, which identifies a peak energy 
iJpk in the ufi^ spectrum, gives reliable fits, though in 5 
cases (class B), the best fit value for _Bpk lies beyond the 
upper limit (150 keV) of the BAT energy band. The un- 
certainties in the fitted parameters are in most cases highly 
correlated among them. 

• The spectra of the remaining 18 bursts (classes C-l-D) 
are well fitted by the simple PL model. For 11 bursts (class 
C) of the 47 GRBs, the best PL-fit value of the photon index 
Q is greater than —2, suggesting that _Epk > 150 keV. The 
other 7 bursts (class D) show a photon index a lower than 
— 2 which implies a decreasing power-law i/fi, spectrum in 
the analysed energy range. In these cases, the peak energy 
could lie in the X-ray range. 

• The fiuence, S.y, and rest-frame isotropic equivalent 
energy, i?T.,iso, as well as their uncertainties, were calculated 
for the 29 (classes A-l-B) GRBs with rehable fits to the CPL 
model. This is currently the largest homogeneous sample of 
long GRBs with determined spectral parameters and E^^iso- 
The mean values of the rest peak energy iSpk ^nd E^^iao in 
the sample are ~ 350 keV and ~ 4 x 10^^ erg, respectively. 

• The i5pk as inferred from the CPL fit and apL as in- 
ferred from the PL fit do correlate. However, this correlation 
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Figure 5. Same as Fig. |4] but for E'^^ and -E-^.iso inferred from spectral fits using the Band model with /3 frozen to —2.3. Note that 
using the Band model the correlation shifts to the higher (lower) E-y isa (^pk) side with respect to the case when the CPL model is 
used. 



is not physical, but is the result of the attempt of the sin- 
gle PL to account for the data points above the peak, with 
smaller flux. As Bpk decreases, a larger fraction of data lies 
above the peak, inducing the single power law to steepen. 
Therefore, the correlation should not be used to infer _Epk 
when knowing only apL- Indeed, we show that i5pk and 
QCPL do not correlate. 

• A correlation between E'p^^ and E-y^i^o (the 'Amati' 
relation) is confirmed for our sample [Eq. ([Sjl]. This correla- 
tion and its scatter are consistent with the ones established 
previously for non-Suii/t bursts, showing that the E^j^-E-y ^iso 
correlation is hardly an artifact of selection effects. The zero- 
point of our correlation (in Log-Log) is larger by ~ O.f dex 
at the sample media n value of j^-y.iBo th an the latest updated 
"Amati" correlation (|Ghirlanda 1120071 ). This difference is ex- 
pected due to the use of a CPL model for describing the 
observed spectra instead of the Band model. 

• When the Band model with the high-energy photon 
index frozen to f3 = —2.3 is fitted to the spectra of our 29 
(classes A+B) GRBs, the zero-point of the resultant E'p^^- 



E-y^iso correlation becomes smaller than in the case of the 
CPL model. For the Band model, the obtained E'p^-E~,^iao 
correlation is given by Eq. ([7|, which is virtually indistinct 
from the 'Amati' correlation established previously for un- 
even observable data-sets from different satellites. 



Although in this work we were able to process ho- 
mogeneously the spectra of a significant fraction of the 
Smft GRBs with known redshift, the sample is still lim- 
ited for conclusive results regarding the physical meaning 
of the correlations among the burst observable properties. 
During the refere eing process of this work, appeared a pa- 
per byQ l|2007bl ). where the author finds evidence of sig- 
nificant change in the 'Amati' correlation parameters with 
re dshift fo r the he terogeneous sample of 48 GRBs compiled 
bv lAmatil (|2006al ). The number of GRBs with known red- 
shifts and full spectral information obtained homogeneously 
should increase in order to attain more conclusive results. 
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APPENDIX A: THE OBSERVED SPECTRA 

In this Appendix section we present the 47 Smft spectra 
(see Table[T} analysed here. As described in §2, the 60 energy 
channels in the range 15-150 keV of the BAT detector are 
used to deconvolve the spectra. The spectral fluxes were 
averaged over the integration time. Tint, determined for each 
burst from its light curve (see §2). In Figs. |XT1 [XU and [Ml 
the CPL model was used to fit the observed spectra. 
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Energy (keV) 

Figure Al. Time-averaged vfi, spectra (red error bars) of the Swift GRBs with known redshift included in our sample. Continuous 
(green) line: best fit curve with the CPL photon model. The spectral flux of each burst was averaged over its iiitegrcition time, T^j^^. 
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Figure A2. Continues Fig. lAll 
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Figure A3. Continues Fig. lAll 



